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The Two-Scale Model and its improved version were used to perform the fit 
to the HERMES data for u (the virtual photon energy) and z (the fraction of v 
| carried by hadron) dependencies of nuclear multiplicity ratios for ir + and ir~ mesons 

I/"") ■ electro-produced on two nuclear targets ( 14 N and 84 Kr). The quantitative criterium 

X 2 was used for the first time to analyse the results of the model fit to the nuclear 
\ multiplicity ratios data. The two-parameter's fit gives satisfactory agreement with 

the HERMES data. Best values of the parameters were then used to calculate the 
v- and z - dependencies of nuclear attenuation for ir°, K + , K~ and p produced on 
I 84 Kr target, and also make a predictions for is, z and the Q 2 (the photon virtuality) 

i ' - dependencies of nuclear attenuation data for those identified hadrons and nuclea, 

gn! that will be published by HERMES. 

& : 

1 Introduction 

Studies of hadron production in deep inelastic semi-inclusive lepton-nucleus scattering 
(SIDIS) offer a possibility to investigate the quark (string, color dipole) propagation in 
dense nuclear matter and the space-time evolution of the hadronization process. It is well- 
known from QCD, that confinement forbids existence of an isolated color charge (quark, 
antiquark, etc.). Consequently, it is clear that after Deep Inelastic Scattering (DIS) of 
lepton on intra-nuclear nucleon, the complicated colorless pre-hadronic system arises. Its 
propagation in the nuclear environment involves processes like multiple interactions with 
the surrounding medium and induced gluon radiation. If the final hadron is formed inside 
the nucleus, the hadron can interact via the relevant hadronic cross section, causing further 
reduction of the hadron yield [1]. QCD at present can not describe the process of quark 
hadronization because of the major role of "soft" interactions. Therefore, the investigation 
of quark hadronization is of basic importance for development of QCD. For this purpose 
we investigate in this paper the Nuclear Attenuation (NA), which is the ratio of the dif- 
ferential multiplicity on nucleus to that on deuterium. At present there exist numerous 
phenomenological models for investigation of the NA problem [2]- [14]. In this work we use 
the Two-Scale Model [4] and its improved version to perform the fit to the HERMES NA 
data [15, 16]. For the fitting purposes we use the more precise part of data, including data 
for v- and z - dependencies of NA of tt + and tt~ mesons on two nuclear targets ( 14 N and 
84 Kr). The v- and z - dependencies of NA for 7r°, K + , K~ and antiproton, produced on 



t 

/ 


\ 


x - 

lc- 
1- 

L - 


- DIIS point 

const, formation length 
yo— yo formation length 

- full hadronization length 












J h 

/ H. 












Ay u 7 /y 








H s 


H,< 

h 4 y 

p< y 


Py 






x x + l c x + l y x + L x 



Figure 1: Space-time structure of hadronization in the string model. The two constituents 
of the hadron are produced at different points. The constituents of the hadron h are created 
at the points P 2 and P 3 . They meet at H 3 to form the hadron. 

84 Kr target we describe with best values of parameters obtained from the above mentioned 
fit. The best set of parameters are used also for prediction of is, z and Q 2 - dependencies 
of NA for the data on those identified hadrons and targets that will be published soon by 
HERMES [17]. The remainder of the paper is organized as follows. In section 2 we briefly 
remind about the Two-Scale Model. In section 3 we discuss the possibility of inclusion 
of the Q 2 -dependence in the Two-Scale Model . In section 4 we describe the scheme we 
used to improve the Two-Scale Model, substituting the step-by-step increase of the string- 
nucleon cross section by a smooth raising function. In section 5 we present the results of 
the model fit to the HERMES data. Our conclusions are given in section 6. 

2 The Two-Scale Model 

The Two-Scale Model is a string model, which was proposed by EMC [4] and used for 
the description of their experimental data. Basic formula is: 

roo roc roo 

R A = 2n bdb dxp(b,x)[l- dx'a str (Ax)p{b,x')] A - 1 (1) 

J J-oo Jx 

where b - impact parameter, x - longitudinal coordinate of the DIS point, x'- longitudi- 
nal coordinate of the string-nucleon interaction point, a str (Ax) - the string-nucleon cross 
section on distance Ax = x'-x from DIS point, p(b, x) - nuclear density function, A - atomic 
mass number. 

The model contains two scale (see Fig. 1): r c (l c ) - constituent formation time (length), 
and Th (lh) - yo-yo formation time (length), 2 yo-yo formation means, that the colorless 

2 in relativistic units (h = c = 1, where h = h/2ir is the Plank reduced constant and c - speed of light) 



system with valence content and quantum numbers of final hadron arises, but without its 
"sea" partons. The simple connection exists between r h and r c 



T h -T c = Zvj K, (2) 

where z= E^/v, Eh and v are energies of final hadron and virtual photon correspon- 
dingly, k - string tension (string constant). Further we will use two different expressions 
for t c . The expression for r c obtained for hadrons containing leading quark [18]: 



t c = {1-z)v/k. (3) 

The expression for average value of r c , which was obtained in [5, 19] in framework of 
the standard Lund model [20]: 

/•OO /'OO 

t c = ldlD c (L,z,l)l dlD c (L,z,l), (4) 

JO JO 

where D C (L, z, 1) is the distribution of the constituent formation length 1 of hadrons 
carrying momentum z. This distribution is: 

D C (L, z, I) = L(l + O + l ° (6(1 -L + zL) + ±±£.)e(l)0{L -zL-l), (5) 

where L = v/n, and parameter C=0.3. The path traveling by string between DIS and 
interaction points is Ax = x'-x. The string-nucleon cross section is: 



a str (Ax) = 6(t c - Ax)a q + 6(r h - Ax)6(Ax - r c )a s + 6(Ax - r h )a h (6) 

where a q , a s and ah are the cross sections for interaction with nucleon of initial string, 
open string (which becomes one of the hadron quarks being looked at) and final hadron 
respectively (see Fig. 2 a)). 

3 Inclusion of the Q 2 -dependence 
in Two-Scale Model. 

The Two-Scale Model [4] does not contain direct Q 2 -dependence and operates with the 
average values of cross sections: 

Oq = °q(Q 2 ); Os = °s(Ql c ), (7) 

where Q 2 is average value of Q 2 obtained in experiment for initial state and Q\ is the 



because after DIS the string radiates gluons and diminishes its virtuality. QCD predicts 
the Q 2 -dependence of string-nucleon cross section in the form [21, 22]: 



a q (Q 2 ) ~ 1/Q 2 ; a s (Q 2 J ~ !/<£■ (8) 
Using this prediction we can write the cross section for initial string as 



° q {Q 2 ) = {Q 2 /Q 2 )° q {Q 2 )- (9) 

In the same way can be written the expression for open string cross section 

= (Q 2 rM> s (Ql), (10) 

where Q 2 c is the virtuality of string for time r c after DIS point, and Q 2 is the same 
for average value of Q 2 . For estimation of ratio Q 2 c /Q 2 c we adopt the scheme given in 
Ref. [23, 24]. In according with this scheme, for the time t the quark decreases its virtuality 
from the initial one, Q 2 , to the value Q 2 (t) 

where z/(t) = v -nt. The calculations shown, that for HERMES kinematics 
(1.2<Q 2 <9.5 GeV 2 and Q 2 =2.5 GeV 2 ), values for ratio Q 2 JQ 2 Tc are close to 1 (for r c in 
form of (3) it changes in region 0.97-^1.04 and in form of (4) in region 0.92^-1.12). This 
means that a s is practically constant. 

4 Improved version of Two-Scale Model 

In the Two-Scale Model the string-nucleon cross section is a function which jumps in points 
Ax = r c and r h . In reality the cross section increases smoothly until it reaches the size 
of hadronic cross section. That is why we need to improve the model in order to obtain 
the smooth increase of the cross section (see Fig. 2). We introduce the parameter c 
(0< c <1) in order to take into account the well known fact, that string starts to interact 
with hadronic cross section soon after creation of the first constituent quark of the final 
hadron, before creation of second constituent. The string-nucleon cross section starts to 
increase from DIS point, and reaches the value of the hadron-nucleon cross section at Ax 
= t. However, in that case one cannot deduce the exact form of a str from perturbative 
QCD, at least in region Ax ~ r. This means, that some model for the shrinkage-expansion 
mechanism has to be invented. We use four versions for a str . Two versions we took from 
Ref. [25] . The first version is based on quantum diffusion: 

a str (Ax) = 0(t - Ax)[a q + (a h - a q )Ax/r] + 6(Ax - r)a h (12) 



where r = r c + cAr, At = - r c , 
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Figure 2: a) The behaviour of the string- nucleon cross section as a function of distance in 
the Two-Scale Model, b) The same as in a) for improved Two-Scale Model with taking 
into account more realistic smoothly increasing string-nucleon cross section. 



a** (Ax) = 6(r - Ax)[a q + (a h - a q )(Ax/rf\ + 6{Ax - r)a h (13) 
We used also two other expressions for a str [2, 6]: 

a** (Ax) = a h - (a h - a q )exp(-Ax/r) (14) 

and: 

a str (Ax) = a h - K - a q )exp(-(Ax/r) 2 ) (15) 

One can easily note that at Ax/t <^ 1 the expressions (14) and (15) turn into (12) and 
(13), correspondingly. 

5 Results 

In this work we have formulated one of possible improvements of the Two-Scale model 
and performed a fit to the HERMES data [15, 16]. Only the data for NA for v- and z - 
dependencies of tt + and n~ mesons on 14 N and 84 Kr nuclea were used for the actual fit. 
Furthermore, NA for v- and z - dependencies of other hadrons produced on 84 Kr target were 
calculated. Also based on the best fit parameters one can make different predictions for u, 



HERMES [17]. The string tension (string constant) was fixed at a static value determined 
by the Regge trajectory slope [24, 26] 



k = 1/(2ttc4) = lGeV/fm (16) 
We use the following Nuclear Density Functions (NDF): 

For 4 He and 14 N we use the Shell Model [27], according to which four nucleons ( two 
protons and two neutrons), fill the s - shell, and other A-4 nucleons are on the p - shell: 

/ 4 2 (A -4) r 2 . . r\ 
P(r) = Po(j + s—A^^'^J' (17) 

where r A —1.31 fm for 4 He and r A =1.67 fm for 14 N. 

For 20 Ne, 84 Kr and 131 Xe we use Woods-Saxon distribution 

p(r) =po/(l + exp((r -r A )/a)). (18) 

The three sets of NDF were used for the fitting with the following corresponding parame- 
ters: 

First set (NDF=1) [28]. 

a = 0.54 fm; r A = (0.978 + 0.0206A 1/3 )A 1/3 fm (19) 
Second set (NDF=2) [29] 

a = 0.54 fm; r A = (1.19A 1/3 - 1.61/A 1/3 ) fm (20) 
Third set (NDF=3) [30] 

a = 0.545 fm; r A = 1.14A 1/3 fm. (21) 
where po are determined from normalization condition: 

1 (22) 



J d 3 rp(r) 



Parameter a is practically the same for all three sets, radius r A for the third set is 
larger approximately on 6% than for second and first sets. From the fit we determined two 
parameters for Two-Scale Model a q and a s . In case of the improved Two-Scale Model the 
fitting parameters are a q and c. 

Determination of the parameter c is represented in Section 4. For fitting we used two 
expressions for r c , which are equations (3) and (4), and five expressions for a str (Ax) (6), 
(12)- (15). The results of the performed fit are presented in Tables 1, 2a and 2b. As we 
have mentioned above only part of HERMES experimental data was used for the fitting 
procedure, including v - and z - dependencies of NA for n + and n~ on 14 N and 84 Kr nuclea. 
For each measured bin the information on the values of z (averaged over the given v bin) 
in case of v dependence, and v in case of z dependence was taken from experimental data. 
Use of this information allows to avoid the problem of additional integration over z and v 
in formulae (1). 




Figure 3: Hadron multiplicity ratio R of charged pions for N and Kr nuclea as a function 
of v (upper panel), z (lower panel). The theoretical curves correspond the calculations in 
Improved Two-Scale Model performed with NDF (17) for 14 N and NDF (19) for 84 Kr and 
a str (12) with r c in form (3) for the values of parameters: <T g =0.46mb, c=0.32. 

In Table 1 the best values for fitted parameters, their errors and x 2 /d.o.f. (N exp =58, 
N par =2. N exp and N par are the numbers of experimental points and fitting parameters 
which were used.) for the Two-Scale Model are represented. Two different expressions for 
r c , and three different sets of parameters for NDF ( 84 Kr) were used. Tables 2a and 2b 
contain the best values for fitted parameters, their errors and x 2 /d.o.f. (N ea , p =58, N par =2) 
for the Improved Two-Scale Model. Four different expressions for a str were used. Only 
difference between Tables 2a and 2b is the form of r c . The results for Two-Scale Model 
(Table 1) are qualitatively close to the results of Ref. [4]. The values of a q <C cr h and a s 
are approximately equal to a h . a q in our case is larger than the same in Ref. [4], because 
Q 2 for HERMES kinematics is smaller than in EMC kinematics. The minimum values for 
X 2 /d.o.f. (best fit) were obtained for the Improved Two-Scale Model with the constituent 
formation time r c in form of (3) (see Table 2a). 

The results for NA, calculated with the best values of fitting parameters for improved 
Two-Scale Model, for v and z dependencies of produced charged pions on 14 N and 84 Kr 
targets are presented on Fig. 3. 



R M h Krypton 



— two par. fit a ) 

- ■ 


-» 






'-. , , , , , , , , , , , , , , , , , , , .«>. 


,,,,,,,,,,,,,,,,, 






'- , , , , , , , , , , , , , , , , , , , .«>. 


^ 






-.. \ .... \ .... \ .... \ .. 


-h) A , 






-.. i .... i .... i .... i .. ^ 


. , . . . , . . . , . . . , Y . 











K + 



K 



5 10 15 20 25 0.2 0.4 0.6 0.8 1.0 

v, GeV z 



Figure 4: Hadron multiplicity ratio R of different species of hadrons produced on 84 Kr 
target [16] as a function of v (left panel) and z (right panel). The curves are calculated 
with the best fit parameters described in the caption of Fig. 3. 



In Fig. 4 one can see the v and z dependencies for all identified hadrons produced on 
84 Kr target. The values of (hadron- nucleon inelastic cross section) used in this work are 
equal to: a n + = cr w - = a n o = Ok- — 20 mb, <Jk+ — 14 mb and ap = 42 mb. The curves 
correspond to the improved Two-Scale model with the best set of parameters. 

In Fig. 5 we present the results of Improved Two-Scale Model in comparison with the 
experimental data for NA of charged hadrons on 63 Cu target [4] performed in region of v 
and Q 2 values higher, than in HERMES kinematics. In order to compare with the EMC 
data we redefined a q to the Q\mc-> according to the expression (9). 

We represent the NA ratio as a function of inverse Q 2 , because of connection of this 
dependence with the Higher Twist effects. Indeed, from the equations (9), (6), (12)-(15) 
and (1) we can conclude, that in first approximation the expansion over the degrees of 
1/Q 2 for NA ratio can be represented in form Ra = a + b/Q 2 , where b is negative. 

One has to note, that for calculation of l/Q 2 -dependence, the cr q (Q 2 ) was used instead 
of a q . Corresponding expression is given by (9). We also take into account nuclear effects 
in deuterium. This means, that instead of a simple formula (1), we use for calculations the 
ratio of (1) for nucleus to the (1) for deuterium. For deuterium as NDF we use Hard Core 




Figure 5: Hadron multiplicity ratio R of charged hadrons for 63 Cu nucleus as a function 
of v ( upper panel) and z (lower panel). The curves are calculated with the best set of 
parameters described in the caption of Fig. 3 
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Figure 6: Hadron multiplicity ratio R of different species of hadrons produced on 4 He target 
as a function of v (left panel), z (central panel) and Q 2 (right panel). The red curves are 
the best fit using improved Two-Scale model with the parameters described in the caption 
of Fig. 3. The black curves correspond to the best fit using the simple Two-Scale model 
with r c defined in (4), NDF (17), a q = 4.2 mb and a s = 16.6 mb 
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Figure 7: The same as described in the caption of the Fig. 5 done for 20 Ne target. 




Figure 8: The same as described in the caption of the Fig. 5 done for 84 Kr target. 
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Figure 9: Descriptive ability of the Two-Scale model and its Improved version: right panel 
- for z, and left panel - for v dependencies of NA. The red lines on both panels correspond 
to the Improved version, the blue ones are for simple Two-Scale model. 

Deuteron Wave Functions from Ref. [31]. 

Using the best set of parameters obtained by fitting the published HERMES data [15, 
16] we calculated the predictions for the new set of the most precise in the world HERMES 
data [17] for 4 He (Fig. 5), 20 Ne (Fig. 6) and 84 Kr (Fig. 7). 

In order to demonstrate the achieved advantages for Imroved Two-Scale model not only 
on the level of obtained x 2 values, one can compare how these two versions are describing 
the NA data for pions on two nuclear targets for z (see right panel of Fig. 8) and v (see 
left panel of Fig. 8) dependencies. It's clearly seen from this plot that being about the 
same for v dependence these two versions remarkable differ for z dependence. 

The last Figure 9 is related to the predictions, done for already presented by HER- 
MES [17] data on 4 He, 20 Ne and 84 Kr targets with the extended kinematics, as well as for 
131 Xe target, on which the data is awaiting soon from the HERMES Collaboration. Two 
set of the best fit parameters were fixed: one marked as a blue curves on Fig. 9 is related to 
the simple Two-Scale Model, next one, marked as a red curves is related to the Improved 
version of the Two-Scale Model. Left panel corresponds to z dependence of NA for pions, 
right panel is related to the v dependence of NA for pions. We can note that again as for 
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Figure 10: Two-Scale model (blue lines) and its improved version (red lines). The predic- 
tions of data on 4 He, 20 Ne, 84 Kr and 131 Xe done for z and v dependencies of NA. 

other nuclear targets, the difference in simple and improved versions is remarkable for Xe 
in z dependence. * 

6 Conclusions. 

• The HERMES data for v- and z - dependencies of nuclear attenuation of 7r + and n~ 
mesons on two nuclear targets ( 14 N and 84 Kr) were used to perform the fit of the 
Two-Scale Model and its Improved Version. 

• Criterion x 2 was use d for the first time to analyse the nuclear attenuation data fit. 

• Two-parameter fit demonstrates satisfactory agreement to the HERMES data. Mini- 
mum x 2 (best fit) was obtained for improved Two-Scale Model, including expressions 
(12) for a str and (3) for r c . The published HERMES data do not give the possibility 
to make a choice between expressions (12)-(15), as well as to prefere definition (3) or 
(4) for r c , because they give close values of x 2 ■ Preferable NDF's are set number one 
and two. 



• More precise data expected from HERMES [17] will provide essentially definite situ- 
ation with the choice of preferable NDF, expressions for a str and r c . 

• In all versions we have obtained a q <C cr h . This indicates that at early stage of 
hadronization process Color Transparency takes place. 



Tc(3) 


r c (4) 


NDF 


a q (mb) 


a s (mb) 


X 2 /d.o.f. 


a q (mb) 


a s (mb) 


X 2 /d.o.f. 


1 


5.3±0.01 


17.1±0.08 


4.3 


4.2±0.01 


16.6±0.07 


2.3 


2 


5.5±0.01 


17.7±0.08 


4.5 


4.3±0.01 


17.3±0.07 


2.4 


3 


5.8±0.010 


18.3±0.08 


4.8 


4.4±0.01 


18.1±0.07 


2.6 



Table 1. The Two-Scale Model. Best values for fitted parameters and x 2 /d.o.f. 

(N ea: p=58, N par =2) 



^(12) 


Cstr(13) 


NDF 


a q (mb) 


c 


X 2 /d.o.f. 


(7 g (mb) 


c 


X 2 /d.o.f. 


1 


0.46±0.02 


0.32±0.03 


1.4 


3.5±0.01 


0.23±0.002 


1.9 


2 


0.62±0.01 


0.31±0.01 


1.7 


3.7±0.01 


0.22±0.02 


2.1 


3 


0.78±0.02 


0.30±0.03 


1.8 


3.9±0.01 


0.21±0.003 


2.3 



<7rtr(14) 




NDF 


cr g (mb) 


c 


X 2 /d.o.f. 


(Tg (mb) 


c 


X 2 /d.o.f. 


1 


l.liO.Ol 


0.15±0.03 


2.1 


3.7±0.01 


0.15±0.02 


2.3 


2 


1.3±0.02 


0.15±0.03 


2.4 


3.9±0.01 


0.14±0.02 


2.6 


3 


1.5±0.02 


0.14±0.03 


2.8 


4.1±0.01 


0.14±0.02 


2.9 



Table 2a. The Improved Two-Scale Model:r c (3). Best values for fitted parameters and 
X 2 /d.o.f. (N exp =58, N par =2). 



^tr(12) 


C S tr( 13) 


NDF 


a q (mb) 


c 


X 2 /d.o.f. 


o q (mb) 


c 


X 2 /d.o.f. 


1 


O.OiO.OOl 


0.56±0.02 


4.6 


0.97±0.01 


0.17±0.002 


1.6 


2 


0.0±0.002 


0.53±0.02 


4.3 


1.0±0.02 


0.17±0.02 


1.5 


3 


0.0±0.002 


0.49±0.006 


4.0 


1.1±0.02 


0.16±0.02 


1.6 



<7rtr(14) 


c str (15) 


NDF 


<jg (mb) 


c 


X 2 /d.o.f. 


(7 9 (mb) 


c 


X 2 /d.o.f. 


1 


0.0±0.001 


0.24±0.02 


3.0 


1.5±0.02 


0.103±0.02 


1.5 


2 


0.0±0.002 


0.21±0.02 


2.9 


1.7±0.02 


0.096±0.02 


1.6 


3 


0.0±0.002 


0.18±0.02 


2.8 


1.8±0.02 


0.089±0.02 


1.8 



Table 2b. The Improved Two-Scale Model:r c (4). Best values for fitted parameters and 
X 2 /d.o.f. (N e;cp =58, N par =2). 



We do not include in consideration NA of protons, because in this case additional 
mechanisms connected with color interaction (string- flip) and final hadron rescattering 
become essential (see for instance Ref. [3, 5]) 

We would like to acknowledge P. Di Nezza and E. Aschenauer for helpful discussions 
and suggestions. We also thank G. Elbakian as well as many other colleagues from the 
HERMES Collaboration for fruitful discussions. 
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The Two-Scale Model and its improved version were used to perform the fit 
to the HERMES data for u (the virtual photon energy) and z (the fraction of v 
| carried by hadron) dependencies of nuclear multiplicity ratios for ir + and ir~ mesons 

I/"") ■ electro-produced on two nuclear targets ( 14 N and 84 Kr). The quantitative criterium 

X 2 was used for the first time to analyse the results of the model fit to the nuclear 
\ multiplicity ratios data. The two-parameter's fit gives satisfactory agreement with 

the HERMES data. Best values of the parameters were then used to calculate the 
v- and z - dependencies of nuclear attenuation for ir°, K + , K~ and p produced on 
I 84 Kr target, and also make a predictions for is, z and the Q 2 (the photon virtuality) 

i ' - dependencies of nuclear attenuation data for those identified hadrons and nuclea, 

gn! that will be published by HERMES. 

& : 

1 Introduction 

Studies of hadron production in deep inelastic semi-inclusive lepton-nucleus scattering 
(SIDIS) offer a possibility to investigate the quark (string, color dipole) propagation in 
dense nuclear matter and the space-time evolution of the hadronization process. It is well- 
known from QCD, that confinement forbids existence of an isolated color charge (quark, 
antiquark, etc.). Consequently, it is clear that after Deep Inelastic Scattering (DIS) of 
lepton on intra-nuclear nucleon, the complicated colorless pre-hadronic system arises. Its 
propagation in the nuclear environment involves processes like multiple interactions with 
the surrounding medium and induced gluon radiation. If the final hadron is formed inside 
the nucleus, the hadron can interact via the relevant hadronic cross section, causing further 
reduction of the hadron yield [1]. QCD at present can not describe the process of quark 
hadronization because of the major role of "soft" interactions. Therefore, the investigation 
of quark hadronization is of basic importance for development of QCD. For this purpose 
we investigate in this paper the Nuclear Attenuation (NA), which is the ratio of the dif- 
ferential multiplicity on nucleus to that on deuterium. At present there exist numerous 
phenomenological models for investigation of the NA problem [2]- [14]. In this work we use 
the Two-Scale Model [4] and its improved version to perform the fit to the HERMES NA 
data [15, 16]. For the fitting purposes we use the more precise part of data, including data 
for v- and z - dependencies of NA of tt + and tt~ mesons on two nuclear targets ( 14 N and 
84 Kr). The v- and z - dependencies of NA for 7r°, K + , K~ and antiproton, produced on 
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Figure 1: Space-time structure of hadronization in the string model. The two constituents 
of the hadron are produced at different points. The constituents of the hadron h are created 
at the points P 2 and P 3 . They meet at H 3 to form the hadron. 

84 Kr target we describe with best values of parameters obtained from the above mentioned 
fit. The best set of parameters are used also for prediction of is, z and Q 2 - dependencies 
of NA for the data on those identified hadrons and targets that will be published soon by 
HERMES [17]. The remainder of the paper is organized as follows. In section 2 we briefly 
remind about the Two-Scale Model. In section 3 we discuss the possibility of inclusion 
of the Q 2 -dependence in the Two-Scale Model . In section 4 we describe the scheme we 
used to improve the Two-Scale Model, substituting the step-by-step increase of the string- 
nucleon cross section by a smooth raising function. In section 5 we present the results of 
the model fit to the HERMES data. Our conclusions are given in section 6. 

2 The Two-Scale Model 

The Two-Scale Model is a string model, which was proposed by EMC [4] and used for 
the description of their experimental data. Basic formula is: 

roo roc roo 

R A = 2n bdb dxp(b,x)[l- dx'a str (Ax)p{b,x')] A - 1 (1) 

J J-oo Jx 

where b - impact parameter, x - longitudinal coordinate of the DIS point, x'- longitudi- 
nal coordinate of the string-nucleon interaction point, a str (Ax) - the string-nucleon cross 
section on distance Ax = x'-x from DIS point, p(b, x) - nuclear density function, A - atomic 
mass number. 

The model contains two scale (see Fig. 1): r c (l c ) - constituent formation time (length), 
and Th (lh) - yo-yo formation time (length), 2 yo-yo formation means, that the colorless 

2 in relativistic units (h = c = 1, where h = h/2ir is the Plank reduced constant and c - speed of light) 



system with valence content and quantum numbers of final hadron arises, but without its 
"sea" partons. The simple connection exists between r h and r c 



T h -T c = Zvj K, (2) 

where z= E^/v, Eh and v are energies of final hadron and virtual photon correspon- 
dingly, k - string tension (string constant). Further we will use two different expressions 
for t c . The expression for r c obtained for hadrons containing leading quark [18]: 



t c = {1-z)v/k. (3) 

The expression for average value of r c , which was obtained in [5, 19] in framework of 
the standard Lund model [20]: 

/•OO /'OO 

t c = ldlD c (L,z,l)l dlD c (L,z,l), (4) 

JO JO 

where D C (L, z, 1) is the distribution of the constituent formation length 1 of hadrons 
carrying momentum z. This distribution is: 

D C (L, z, I) = L(l + O + l ° (6(1 -L + zL) + ±±£.)e(l)0{L -zL-l), (5) 

where L = v/n, and parameter C=0.3. The path traveling by string between DIS and 
interaction points is Ax = x'-x. The string-nucleon cross section is: 



a str (Ax) = 6(t c - Ax)a q + 6(r h - Ax)6(Ax - r c )a s + 6(Ax - r h )a h (6) 

where a q , a s and ah are the cross sections for interaction with nucleon of initial string, 
open string (which becomes one of the hadron quarks being looked at) and final hadron 
respectively (see Fig. 2 a)). 

3 Inclusion of the Q 2 -dependence 
in Two-Scale Model. 

The Two-Scale Model [4] does not contain direct Q 2 -dependence and operates with the 
average values of cross sections: 

Oq = °q(Q 2 ); Os = °s(Ql c ), (7) 

where Q 2 is average value of Q 2 obtained in experiment for initial state and Q\ is the 



because after DIS the string radiates gluons and diminishes its virtuality. QCD predicts 
the Q 2 -dependence of string-nucleon cross section in the form [21, 22]: 



a q (Q 2 ) ~ 1/Q 2 ; a s (Q 2 J ~ !/<£■ (8) 
Using this prediction we can write the cross section for initial string as 



° q {Q 2 ) = {Q 2 /Q 2 )° q {Q 2 )- (9) 

In the same way can be written the expression for open string cross section 

= (Q 2 rM> s (Ql), (10) 

where Q 2 c is the virtuality of string for time r c after DIS point, and Q 2 is the same 
for average value of Q 2 . For estimation of ratio Q 2 c /Q 2 c we adopt the scheme given in 
Ref. [23, 24]. In according with this scheme, for the time t the quark decreases its virtuality 
from the initial one, Q 2 , to the value Q 2 (t) 

where z/(t) = v -nt. The calculations shown, that for HERMES kinematics 
(1.2<Q 2 <9.5 GeV 2 and Q 2 =2.5 GeV 2 ), values for ratio Q 2 JQ 2 Tc are close to 1 (for r c in 
form of (3) it changes in region 0.97-^1.04 and in form of (4) in region 0.92^-1.12). This 
means that a s is practically constant. 

4 Improved version of Two-Scale Model 

In the Two-Scale Model the string-nucleon cross section is a function which jumps in points 
Ax = r c and r h . In reality the cross section increases smoothly until it reaches the size 
of hadronic cross section. That is why we need to improve the model in order to obtain 
the smooth increase of the cross section (see Fig. 2). We introduce the parameter c 
(0< c <1) in order to take into account the well known fact, that string starts to interact 
with hadronic cross section soon after creation of the first constituent quark of the final 
hadron, before creation of second constituent. The string-nucleon cross section starts to 
increase from DIS point, and reaches the value of the hadron-nucleon cross section at Ax 
= t. However, in that case one cannot deduce the exact form of a str from perturbative 
QCD, at least in region Ax ~ r. This means, that some model for the shrinkage-expansion 
mechanism has to be invented. We use four versions for a str . Two versions we took from 
Ref. [25] . The first version is based on quantum diffusion: 

a str (Ax) = 0(t - Ax)[a q + (a h - a q )Ax/r] + 6(Ax - r)a h (12) 



where r = r c + cAr, At = - r c , 
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Figure 2: a) The behaviour of the string- nucleon cross section as a function of distance in 
the Two-Scale Model, b) The same as in a) for improved Two-Scale Model with taking 
into account more realistic smoothly increasing string-nucleon cross section. 



a** (Ax) = 6(r - Ax)[a q + (a h - a q )(Ax/rf\ + 6{Ax - r)a h (13) 
We used also two other expressions for a str [2, 6]: 

a** (Ax) = a h - (a h - a q )exp(-Ax/r) (14) 

and: 

a str (Ax) = a h - K - a q )exp(-(Ax/r) 2 ) (15) 

One can easily note that at Ax/t <^ 1 the expressions (14) and (15) turn into (12) and 
(13), correspondingly. 

5 Results 

In this work we have formulated one of possible improvements of the Two-Scale model 
and performed a fit to the HERMES data [15, 16]. Only the data for NA for v- and z - 
dependencies of tt + and n~ mesons on 14 N and 84 Kr nuclea were used for the actual fit. 
Furthermore, NA for v- and z - dependencies of other hadrons produced on 84 Kr target were 
calculated. Also based on the best fit parameters one can make different predictions for u, 



HERMES [17]. The string tension (string constant) was fixed at a static value determined 
by the Regge trajectory slope [24, 26] 



k = 1/(2ttc4) = lGeV/fm (16) 
We use the following Nuclear Density Functions (NDF): 

For 4 He and 14 N we use the Shell Model [27], according to which four nucleons ( two 
protons and two neutrons), fill the s - shell, and other A-4 nucleons are on the p - shell: 

/ 4 2 (A -4) r 2 . . r\ 
P(r) = Po(j + s—A^^'^J' (17) 

where r A —1.31 fm for 4 He and r A =1.67 fm for 14 N. 

For 20 Ne, 84 Kr and 131 Xe we use Woods-Saxon distribution 

p(r) =po/(l + exp((r -r A )/a)). (18) 

The three sets of NDF were used for the fitting with the following corresponding parame- 
ters: 

First set (NDF=1) [28]. 

a = 0.54 fm; r A = (0.978 + 0.0206A 1/3 )A 1/3 fm (19) 
Second set (NDF=2) [29] 

a = 0.54 fm; r A = (1.19A 1/3 - 1.61/A 1/3 ) fm (20) 
Third set (NDF=3) [30] 

a = 0.545 fm; r A = 1.14A 1/3 fm. (21) 
where po are determined from normalization condition: 

1 (22) 



J d 3 rp(r) 



Parameter a is practically the same for all three sets, radius r A for the third set is 
larger approximately on 6% than for second and first sets. From the fit we determined two 
parameters for Two-Scale Model a q and a s . In case of the improved Two-Scale Model the 
fitting parameters are a q and c. 

Determination of the parameter c is represented in Section 4. For fitting we used two 
expressions for r c , which are equations (3) and (4), and five expressions for a str (Ax) (6), 
(12)- (15). The results of the performed fit are presented in Tables 1, 2a and 2b. As we 
have mentioned above only part of HERMES experimental data was used for the fitting 
procedure, including v - and z - dependencies of NA for n + and n~ on 14 N and 84 Kr nuclea. 
For each measured bin the information on the values of z (averaged over the given v bin) 
in case of v dependence, and v in case of z dependence was taken from experimental data. 
Use of this information allows to avoid the problem of additional integration over z and v 
in formulae (1). 




Figure 3: Hadron multiplicity ratio R of charged pions for N and Kr nuclea as a function 
of v (upper panel), z (lower panel). The theoretical curves correspond the calculations in 
Improved Two-Scale Model performed with NDF (17) for 14 N and NDF (19) for 84 Kr and 
a str (12) with r c in form (3) for the values of parameters: <T g =0.46mb, c=0.32. 

In Table 1 the best values for fitted parameters, their errors and x 2 /d.o.f. (N exp =58, 
N par =2. N exp and N par are the numbers of experimental points and fitting parameters 
which were used.) for the Two-Scale Model are represented. Two different expressions for 
r c , and three different sets of parameters for NDF ( 84 Kr) were used. Tables 2a and 2b 
contain the best values for fitted parameters, their errors and x 2 /d.o.f. (N ea , p =58, N par =2) 
for the Improved Two-Scale Model. Four different expressions for a str were used. Only 
difference between Tables 2a and 2b is the form of r c . The results for Two-Scale Model 
(Table 1) are qualitatively close to the results of Ref. [4]. The values of a q <C cr h and a s 
are approximately equal to a h . a q in our case is larger than the same in Ref. [4], because 
Q 2 for HERMES kinematics is smaller than in EMC kinematics. The minimum values for 
X 2 /d.o.f. (best fit) were obtained for the Improved Two-Scale Model with the constituent 
formation time r c in form of (3) (see Table 2a). 

The results for NA, calculated with the best values of fitting parameters for improved 
Two-Scale Model, for v and z dependencies of produced charged pions on 14 N and 84 Kr 
targets are presented on Fig. 3. 
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Figure 4: Hadron multiplicity ratio R of different species of hadrons produced on 84 Kr 
target [16] as a function of v (left panel) and z (right panel). The curves are calculated 
with the best fit parameters described in the caption of Fig. 3. 



In Fig. 4 one can see the v and z dependencies for all identified hadrons produced on 
84 Kr target. The values of (hadron- nucleon inelastic cross section) used in this work are 
equal to: a n + = cr w - = a n o = Ok- — 20 mb, <Jk+ — 14 mb and ap = 42 mb. The curves 
correspond to the improved Two-Scale model with the best set of parameters. 

In Fig. 5 we present the results of the Two-Scale Model and its improved version 
in comparison with the experimental data for NA of charged hadrons on 63 Cu target [4] 
performed in region of v and Q 2 values higher, than in HERMES kinematics. In order to 
compare with the EMC data we redefined a q to the Qemci according to the expression 
(9). 

We represent the NA ratio as a function of inverse Q 2 , because of connection of this 
dependence with the Higher Twist effects. Indeed, from the equations (9), (6), (12)-(15) 
and (1) we can conclude, that in first approximation the expansion over the degrees of 
1/Q 2 for NA ratio can be represented in form Ra = a + b/Q 2 , where b is negative. 

One has to note, that for calculation of l/Q 2 -dependence, the <T g (Q 2 ) was used instead 
of o q . Corresponding expression is given by (9). We also take into account nuclear effects 
in deuterium. This means, that instead of a simple formula (1), we use for calculations the 




Figure 5: Hadron multiplicity ratio R of charged hadrons for Cu nucleus as a function 
of v (upper panel) and z (lower panel). The dashed curves correspond to the Two-Scale 
Model, the solid ones to the improved version. The solid curves are calculated with the 
best set of parameters described in the caption of Fig. 3. The dashed curves are calculated 
with r c (4), NDF (19), a q = 4.2 mb nad a s = 16.6 mb. 
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Figure 6: Hadron multiplicity ratio R of different species of hadrons produced on 4 He target 
as a function of v (left panel), z (central panel) and Q 2 (right panel). The solid curves are 
the best fit using improved Two-Scale model with the parameters described in the caption 
of Fig. 3. The dashed curves correspond to the best fit using the simple Two-Scale model 
with r c defined in (4), NDF (17), a q = 4.2 mb and a s = 16.6 mb 
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Figure 7: The same as described in the caption of the Fig. 6 done for 20 Ne target. 




Figure 8: The same as described in the caption of the Fig. 6 done for 84 Kr target. 
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Figure 9: Descriptive ability of the Two-Scale model and its Improved version: right panel 
- for z, and left panel - for v dependencies of NA. The solid lines on both panels correspond 
to the Improved version, the dashed ones are for simple Two-Scale model. 

ratio of (1) for nucleus to the (1) for deuterium. For deuterium as NDF we use Hard Core 
Deuteron Wave Functions from Ref. [31]. 

Using the best set of parameters obtained by fitting the published HERMES data [15, 
16] we calculated the predictions for the new set of the most precise in the world HERMES 
data [17] for 4 He (Fig. 6), 20 Ne (Fig. 7) and 84 Kr (Fig. 8). 

In order to demonstrate the achieved advantages for Imroved Two-Scale model not only 
on the level of obtained \ 2 values, one can compare how these two versions are describing 
the NA data for pions on two nuclear targets for z (see right panel of Fig. 9) and v (see 
left panel of Fig. 9) dependencies. It's clearly seen from this plot that being about the 
same for v dependence these two versions remarkable differ for z dependence. 

The last Figure 10 is related to the predictions, done for already presented by HER- 
MES [17] data on 4 He, 20 Ne and 84 Kr targets with the extended kinematics, as well as for 
131 Xe target, on which the data is awaiting soon from the HERMES Collaboration. Two 
set of the best fit parameters were fixed: one marked as a dashed curves on Fig. 10 is 
related to the simple Two-Scale Model, next one, marked as a solid curves is related to 
the Improved version of the Two-Scale Model. Left panel corresponds to z dependence of 
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Figure 10: Two-Scale model (dashed lines) and its improved version (solid lines). The 
predictions for data on 4 He, 20 Ne, 84 Kr and 131 Xe done for z and v dependencies of NA. 

NA for pions, right panel is related to the v dependence of NA for pions. We can note 
that again as for other nuclear targets, the difference in simple and improved versions is 
remarkable for Xe in z dependence. 

6 Conclusions. 

• The HERMES data for v- and z - dependencies of nuclear attenuation of ir + and 7r~ 
mesons on two nuclear targets ( 14 N and 84 Kr) were used to perform the fit of the 
Two-Scale Model and its Improved Version. 

• Criterion x 2 was used for the first time to analyse the nuclear attenuation data fit. 

• Two-parameter fit demonstrates satisfactory agreement to the HERMES data. Mini- 
mum x 2 (best fit) was obtained for improved Two-Scale Model, including expressions 
(12) for a str and (3) for r c . The published HERMES data do not give the possibility 
to make a choice between expressions (12)-(15), as well as to prefere definition (3) or 



(4) for r c , because they give close values of x 2 - Preferable NDF's are set number one 
and two. 

• More precise data expected from HERMES [17] will provide essentially definite situ- 
ation with the choice of preferable NDF, expressions for a str and r c . 

• In all versions we have obtained a q <C cr h . This indicates that at early stage of 
hadronization process Color Transparency takes place. 



Tc(3) 


r c (4) 


NDF 


a q (mb) 


<r s (mb) 


X 2 /d.o.f. 


a q (mb) 


a s (mb) 


X 2 /d.o.f. 


1 


5.3±0.01 


17.1±0.08 


4.3 


4.2±0.01 


16.6±0.07 


2.3 


2 


5.5±0.01 


17.7±0.08 


4.5 


4.3±0.01 


17.3±0.07 


2.4 


3 


5.8±0.01 


18.3±0.08 


4.8 


4.4±0.01 


18.1±0.07 


2.6 



Table 1. The Two-Scale Model. Best values for fitted parameters and x 2 /d.o.f. 

(N ea: p=58, N par =2) 





ow(13) 


NDF 


cr q (mb) 


c 


X 2 /d.o.f. 


(T q (mb) 


c 


X 2 /d.o.f. 


1 


0.46±0.02 


0.32±0.03 


1.4 


3.5±0.01 


0.23±0.002 


1.9 


2 


0.62±0.01 


0.31±0.01 


1.7 


3.7±0.01 


0.22±0.02 


2.1 


3 


0.78±0.02 


0.30±0.03 


1.8 


3.9±0.01 


0.21±0.003 


2.3 



<7rtr(14) 


ow(15) 


NDF 


o q (mb) 


c 


X 2 /d.o.f. 


o- 9 (mb) 


c 


X 2 /d.o.f. 


1 


l.liO.Ol 


0.15±0.03 


2.1 


3.7±0.01 


0.15±0.02 


2.3 


2 


1.3±0.02 


0.15±0.03 


2.4 


3.9±0.01 


0.14±0.02 


2.6 


3 


1.5±0.02 


0.14±0.03 


2.8 


4.1±0.01 


0.14±0.02 


2.9 



Table 2a. The Improved Two-Scale Model:r c (3). Best values for fitted parameters and 
X 2 /d.o.f. (N exp =58, N par =2). 



^tr(12) 


C S tr( 13) 


NDF 


a q (mb) 


c 


X 2 /d.o.f. 


o q (mb) 


c 


X 2 /d.o.f. 


1 


O.OiO.OOl 


0.56±0.02 


4.6 


0.97±0.01 


0.17±0.002 


1.6 


2 


0.0±0.002 


0.53±0.02 


4.3 


1.0±0.02 


0.17±0.02 


1.5 


3 


0.0±0.002 


0.49±0.006 


4.0 


1.1±0.02 


0.16±0.02 


1.6 



<7rtr(14) 


c str (15) 


NDF 


<jg (mb) 


c 


X 2 /d.o.f. 


(7 9 (mb) 


c 


X 2 /d.o.f. 


1 


0.0±0.001 


0.24±0.02 


3.0 


1.5±0.02 


0.103±0.02 


1.5 


2 


0.0±0.002 


0.21±0.02 


2.9 


1.7±0.02 


0.096±0.02 


1.6 


3 


0.0±0.002 


0.18±0.02 


2.8 


1.8±0.02 


0.089±0.02 


1.8 



Table 2b. The Improved Two-Scale Model:r c (4). Best values for fitted parameters and 
X 2 /d.o.f. (N e;cp =58, N par =2). 



We do not include in consideration NA of protons, because in this case additional 
mechanisms connected with color interaction (string- flip) and final hadron rescattering 
become essential (see for instance Ref. [3, 5]) 

We would like to acknowledge P. Di Nezza and E. Aschenauer for helpful discussions 
and suggestions. We also thank G. Elbakian as well as many other colleagues from the 
HERMES Collaboration for fruitful discussions. 
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